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Abstract 

A multiparticle system produced by a large number of independent sources is described 
by a gaussian density matrix W. All theoretical approaches to Bose-Einstein Correlations 
Cn in high energy physics use this form for W. One of the most salient consequences of this 
form is the fact that all higher order (n > 2) moments of the current distribution can be 
expressed in terms of the first two. We test this property by comparing the data on C 2 (Q^), 
and C 4 (Q^) from 7r+p and A+p reactions at 250 GeVjc with the predictions of 
a general quantum statistical space-time approach. Even a simplified version of such an 
approach involving only 4 (instead of the total of 10) independent parameters (proper-time, 
correlation-length, transverse radius of chaotic source, chaoticity) can account for the data. 
Previous attempts along these lines, which did not use the space-time approach, met with 
difficulties. 
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Bose-Einstein Correlations (BEC) are the basis of an experimental method for the determina¬ 
tion of sizes and lifetimes of sources in particle and nuclear physics. This knowledge is essential 
for an understanding of the dynamics of strong interactions. 

A particularly important aspect of BEC is represented by higher order correlations because of 
the prediction, which follows from the gaussian density matrix assumption, that all higher order 
moments of the current distribution can be described in terms of the first two. The gaussian form 
(in the coherent state representation) of the density matrix Q is a fundamental assumption of 
BEC and follows from the central limit theorem for a large number of independent sources, which 
are expected to act in a high energy reaction. Furthermore, higher order correlations provide 
important constraints on the space-time form of the sources, their dynamics (expansion) and 
chaoticity. 

Experimentally correlations of three and more particles have been studied in the last years in 
[||-|^] and more recently attempts have been made to analyze these correlations in terms of sim¬ 


plified models Il^]-|1T| without clear space-time implications for the emitting source (sources). 
Usually gaussian or exponential forms for the correlation functions in momentum space are 
postulated. In contrast to this, the approach used in this work starts with the space-time char¬ 
acteristics of an expanding source and the space-time form of the correlators within the classical 
current formalism. The dependence of Cn on the four-momentum difference {Q) follows after 


explicit integration over space-time variables |16|. 

The aim of this investigation is to use the higher order correlations NA22 data to test the 
validity of the gaussian density matrix assumption, within the space-time approach to BEC. 


^The mathematical form of the density matrix must not be confused with the form of the correlator or of the 
space-time distribution of the source, cf.below. 
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1 The General Formalism 


In quantum mechanics, a multiparticle production process is described in terms of a density 
matrix W which characterizes the final state of the system. From the density matrix, all n- 
particle distributions can be determined, and conversely, a measurement of these distributions 
yields information about the density matrix of the multiparticle system. The n-particle inclusive 
distribution is defined through the creation and annihilation operator a|(k) and Oj(k) of a 
particle of momentum k (i labels internal degrees of freedom): 


pj,i-*"(ki,...,k„) = 
1 




where. 


diOi = 


d^kj 


is the invariant volume in momentum space. 

The general n-particle correlation function is defined as 


rh-ir.(X, W ^ ^ .. C^l-^”(ki,...,kn) 

" ^ piHkl)...pl"(kn) pr(kl)...pr(k„) ■ 


( 1 ) 

( 2 ) 

( 3 ) 


In order to determine the density matrix for a given reaction from first principles, one would 
have to specify the initial state of the projectile and target and then apply the S matrix to this 
state. In general this is not possible. One way to proceed is to parametrize W according to a 
reasonable phenomenological description of the system. For this, one uses the external source 
(current) formalism. In this approach particle sources are treated as external classical currents, 
and their fluctuations are described by a gaussian distribution. This last choice can be justified 
by the fact that, if one has a superposition of N independent sources, the gaussian form follows 
from the central limit theorem in the limit of large N. 

In the following we will also consider correlation functions as functions of ^ 

^ As usual, Oi = J2i<j lij’ lij = n>2 and i,j = 1,..., n. 
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with (i, j = 1, 


(J (Q^) — I I ) 


( 4 ) 


In{Q^) = I dWi... I dwn cy^-{ku -, K) 6[Q^ + J2 - %)'] 


(5) 


ij=l(i<i) 


lUQ^) = / (h - kj] 

i,j = l{i<j) 


( 6 ) 


Here we give a brief summary of the derivation of Bose-Einstein correlation functions in the 
current formalism ||T7| . 

The current can in general be written as the sum of a chaotic and a coherent component, 
J{x) = Jchaotic{x) + Jcoherent{x) ■ The Gaussian Current distribution is completely specified by 
its hrst two moments: I{x) =< J{x) >=Jcoherent{x) and the two-current correlator 
D{x,y) =< J{x)J{y) > - < J{x) X J{y) >. 

I{x) and D{x,y) can be parametrized as I{x) oc fdx) and D{x,y) oc fch{x)C{x — y)fch{y), 
where fdx) and fchix) are the space-time distributions of the coherent and the chaotic compo¬ 
nents of the source. The primordial correlator C{x — y) reflects intrinsic dynamical properties 
of the source. It contains some characteristic length (or time) scales L, so-called correlation 
lengths (for a system in thermal equilibrium the correlation length can be related to the inverse 
of the temperature). 

In the general case of a partially coherent source, the single inclusive distributions of pions 
can be expressed also as a sum of a chaotic and coherent component (i = ,0 denotes the 

charge): 


1 da^ 
a dio 


a duj 


chaotic 


1 da^ 
a dio 


coherent 


(7) 


where, 


a doj 


chaotic 


D{k), 


( 8 ) 
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1 da^ 
a doj 


= \im^- 


(9) 


I coherent 

I{k) and D{k) are the on-shell Fourier transforms of I{x) and D{x,y), respectively. 
In general, the chaoticity parameter will be momentum-dependent: 

D{k) 


p{k) = 


D{k) + \I{k)\^' 


( 10 ) 


To write down the correlations functions in a concise form one introduces the normalized 
current correlator: 

, _ D{kr, kg') 

[D{kr,kr).D{ks,ks)V/^ 

where the indices r,s label the particles. Since d{kr,ks) is in general a complex number, one 
may prefer to express the correlation functions in terms of the magnitudes and the phases of 

drs ) 



T 

rs 


(12) 

^rs 

= (lf^{kr,ks) = Q.Tgd{kr,ks) 

(13) 


and the phase of the coherent component 

cPl° = r{kr)=S.TgI{kr). (14) 

With this notation we write the chaoticity parameter: pr =p{kr). 

The normalized cumulant correlation functions Hn for identical charged particles are: 

H:^^{ki,k2) = 2y^pi(l - pl)p2{l - P2)Ti2Cos{4r^ - -t- (/>“) + P 1 P 2 T 12 (15) 

H+++{ki,k2, fcs) = lpiP2P3TuT23niCOsicPf2 + <^>^3 + 

~k\l Pl(l “ pl)p24'3(l “ P3)Tl2T23COs{(j)i2 + 4>23 + 4'T ~ 

-|-permutation of 1,2,3 (16) 
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Hl^~^^{ki,k 2 ,h,k 4 ) = 


-^PlP2P3PiTi2T23T34T4iCOs{(p'l2 + </>23 + <^'34 + </'4l) 

+ \/pi(l ~ P^)P2P3P‘ii^ — P4:)Tl2T23T34COs{(pl2 + ^23 + ^3A + 4>T ~ 4>T) 

+permutation of 1,2,3,4 


( 17 ) 


The correlation functions Cn (3) can be expressed in terms of by using the relation 


between the corresponding generating functionals |17]: 




(18) 


Cl^'^^'^%ki,k2,k3) = 

1 + fcs) + 77*2’'^(^2, h) + Hl^’^^iks, ki) + Hi^’^^’^^ki,k2, k^) 


(19) 


^*1,*2,*3,^4 


( ki,k2,k3,k4) = 

1 + 77^'^(A:i,fe2) + 77^^’'^(A;i,A:3) + 77‘^’'"(A:i,A:4) + 77?’*^(A;2,A:3) 

+ Hi^'^^{k2,k4) + Hl^''^{k3,k4) 

+ 77^'^(A:i, fe2).77?’'^(fe2, h) + Hl^^^\k4M)-HT'\k2, ^ 4 ) 

+ Hi^^^\k4,k4).Hl^^'\k2,k3) 

+ Wi'^^'^^{kuk2, h) + fc2, k4) + 77*i’*^’*^(A:i, fcg, ^ 4 ) 


( 20 ) 


As one can see, all correlation functions depend only on the functions T^^s, 4>r^s 4’r° which 

will be specified in the next section for an expanding source. 

For the space-time description of an expanding source is useful to define variables r, r] and 


xii: 


r = 



- XII 



Xo + X 
Xo — X 


( 21 ) 
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where r is the proper-time, xy the coordinate in the direction of the collision axis and r] the 
space-time rapidity. Here we will consider invariance under boosts of the coordinate frame in the 
longitudinal direction, which corresponds to the assumption that the single inclusive distribution 
in rapidity is flat. 

The space-time distribution of the chaotic and coherent source, as well as the primordial cor¬ 
relator 0, are expressed in terms of 10 parameters: ro,ch and ro,co are the proper time coordinates 
of the chaotic and the coherent source; 6Tch and 6tco their widths in proper-time; Rch and Rco 
are the transverse radii; L±_ and are the correlation lengths and is the correlation time. 
The tenth one is the chaoticity parameter. 

To reduce the number of parameters we assume now that the widths in proper-time of both 
sources are vanishing {6Tch = Stco = 0); with this choice the results do not depend on the 
correlation time L^-. The space-time distributions of the chaotic and coherent sources are then 
parametrized as: 


fch ~ 5(r - ro,c/i)exp (22) 

fco ~ S{t - To,co)exp 

Note that the rj dependence of fi{x) was neglected in equations (22) and (23). This corre¬ 
sponds to a boost-invariant ansatz of the source expansion. 

The model contains now 7 independent parameters: tq^c/d "To.co, Rch-, Rco-, T_l, Lr ^ and the 
chaoticity parameter pq 2. 

®The parametrization for the primordial correlator, in the case of an expanding source, has to take into account 
that each source element is characterized not only by a correlation length but also by a four-velocity. Effects of the 
geometry of the source are considered by introducing the space-time distributions of the chaotic and the coherent 
component, fch{x) and fco{x) respectively. 

^The relative contributions of the chaotic and coherent component are determined by fixing the value of the 
(momentum-dependent) chaoticity parameter p at same arbitrary scale, e.g., po =p(k = 0) (cf. Eq. (32)). 
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We will define, for convenience: 


b = 


' 0,c/i 


p2 T 2 


Rlk + ri 


and 712 = 


T0,cfe(mi_L - m2±) 

L2mi_Lm2± 


The single inclusive distribution will be a sum of a chaotic and a coherent term: 


cj 

Sch{k) 


Scoik^ 


\poSch{k) + (1 -po)sco{k)]E 


1 d?a 
a d?k 


A:=0 


’’’"K 

- exp 

m± 



rriT, 

- exp 

m_L 


klRlo \^ 


(24) 


(25) 

(26) 
(27) 


where m± is the transverse mass of the pions emitted and the momentum dependence of chaotic- 
ity parameter takes the form (r, s = 1, 2, 3,4): 


Pr = p{kr) = ^ , Ar =Po + {l-po)Sr 

Xi - Y ' 


Srs = exp 


{kj^ + k-^^^ml-Rl) 


The magnitudes and phases of Trs, 4>rs s-iid 4’?’ are: 

b{yr-ysf {kr±-k‘l^){Rl^- 


Trs = (1 + 7 rJ ^^'^exp 


1 + 7, 


(28) 

(29) 


(30) 


ich 

Ks 


h^r 


1 +7r^ 


■{yr-Vsf - TQ^ch{mrX “ rUax) - - arctan 7 ^^ 


(31) 


(K° = -T0,comr± 


(32) 


We proceed now to the the calculation of the correlation functions. 


2 Calculations and Comparison with Experimental Data 


The procedure for calculations consists in: 

a) The phase-space is generated by Monte Carlo routines which take into account the exper¬ 
imental detection conditions mentioned in (defined by \xf\ < 0.5) and references 
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quoted there |j. 

b) The simulated produced pious are registered with a determined y, k± and azimuth angle 4> 
and from them one calculates the correspondent value. The simulated particles are then 
selected in bins, as described in the experimental papers quoted above. 

c) The calculation of the functions Trs, (l>'^ and defined for an expanding source model is 
performed. From these results one is able to calculate the correlation function (and its integral 
for each defined bin), as determined by equations (4)-(6) and (18)-(20). 

d) We use this procedure for second, third and fourth order correlation calculations, C' 2 (< 5 ^), C 3 (Q^) 
and 6*4 (Q^). 

To simplify the calculations we reduce the number of free parameters using the criteria 


presented in |16|. 

We assume at first that the chaotic and coherent components have the same transverse 
momentum spectrum, i.e., we choose Rco = Rl (Ifpo = 1 ) Rco is obviously not present anymore 
in the formalism.). Rl is constrained by the relation | Rl = {^/-k/ 2)< kt We also take 


Tchfi — 'TcOjO — TQ- 

We are thus left with only four parameters: tq, Rch, L-q and po- 

These parameters have been investigated in the intervals]]: 0.1 < po < 1.0; 0.5 < Rch < 3.5; 
1.0 < To < 2.5 and 0.1 < < 1.0. 

We will use the x^/ndf {ndf: number of degrees of freedom) method to evaluate the quality of 
the fits. We are looking for a common fit for all C'j(Q^) using one set of parameters {po,Rch,To,Lri). 
The procedure to select this one set for an overall fit was: 

®We have also investigated the sensitivity of the fits to the phase space varying the acceptance conditions. 
For this purpose simulations for a different experimental window —2 < y < 2, 0.125G'eF < kt < 1.5GeV were 
performed. We found that in the range of < O.SGeF^ the third and fourth order correlation functions are 
more sensitive to the detections conditions than the second order correlation function. 

®< fct > is taken from experiment. 

^The steps used in the calculation to adjust the parameters are ARch = Apo = Aro = ALn = 0.1. 


9 





1) Calculation of for each Ci{Q‘^) {i = 2,3,4) function for all parameter sets to delimit 
sectors in the parameters phase-space to be investigated. The groups of parameters which gave 
X^/ndf> 3 for any calculated function were eliminated. 

2) Search for regions of intersection in the parameters space for an overall description of the 
data. These regions were defined by minimizing /"^df for all three functions simultaneously 
and imposing x?/nd/< 1-5- 

From 1) we concluded that there is a large number of parameter combinations which permit 
acceptable fits for 6*2 (Q^). The situation is quite different when one extends the search to 
C' 3 (Q 2 ) and C 4 (Q 2 ). 

Table 1 shows the groups of parameters we would like to comment on primarily. In all these 
groups po = 1- This table also contains the /^df for each calculation. Figure 1 shows the 
results for the best overall fit ^ corresponding to the parameters set N2 P^. 

Some results for po ^ 0 should be commented. Calculations using pQ = 0.3 and Rch = 1.7/m 
can produce good fits for C' 2 (Q^) (xi/^^f ^ 1-3) fo < Ln < 0.7/m/c and in the entire 

To interval, but the minimum values for Ca and (74 are around 2.0 in the (very small) 

region of intersection of parameters. The behaviour of the C 2 function does not change even if 
the radius varies up to 2.5 fm, but then there does not exist an intersection region of parameters 
(with xl/ndf< 1.5) for Cj, and C^. The calculations using po = 0.7 presented almost the same 
behaviour. 

We found a clear tendency towards a common value of = 0.3/m/c. The prefered range 

of the radius is 1.7fm < R^h < 2.2/m. The range of tq seems to be connected with that of the 
®The best fits do not necessarily present the lowest value for x^/nd/ of each function, since we look for regions 

in the parameter phase-space which minimize the x^/ud/ for all three functions. 

®We normalized the calculated correlation functions by choosing Ci{Q^) = 1 for = 2.0GeF^ {i = 2,3,4). 

The normalization factor is very close to one. 

^°For completeness, calculations have been done also for Cn{Q) besides Cn{Q^). The Cn{Q) functions did not 
lead to new conclusions about the choice of parameters, they were in general agreement with the behaviour of 

G„(g"). 
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radius. One gets the best fits for R^h = 1.7/m with tq = l.lfm/c and for Rch = 2.2/m with 
To = l.bfmic (see Table 1) 

As already mentioned, the number of parameter sets which provide an acceptable fit for 
C 2 {Q'^) is much bigger than the number of parameter sets which provide an acceptable general 
fit. This is exemplified in Figure 2. This figure should also also clarify the method used in this 
work, which differs from the one used in ||2^, where the parameters obtained from a fit of C 2 
data were used to “predict” C^. 

For the results in Figure 2 we used: po = 0.7, Rch = 1.8/m, tq = 2.2/m/c and = 0.5fm/c. 
The corresponding /ndf values for C 2 , and functions are: 1.30, 7.05, 1.83. 


3 Discussion and Conclusions 

Our main conclusion is that in the frame of a model based on quantum statistical principles and a 
space-time picture of an expanding source, a general description of the higher order correlations 
data with the same parameters as those appearing in the first two correlation orders, is 
possible. The assumption of a gaussian form for the density matrix is consistent with the data. 

A similar but weaker conclusion was reached in |]l^ . In |T^ it had been shown that the 
gaussian form of the density matrix is robust enough to resist attempts of falsification § 

On the other hand in and [^] the momentum-space approach was used, the currents 
We have also calculated the statistical errors associated with our fits. This was done by fixing, at a given 
run, 3 parameters out of the total 4 and searching until an increase of x? value by one unity was obtained. The 
results for the set N2 (poRchRv^'^o) are: (0.6 - 1.0, 2.2, 0.3,1.5); (1.0,1.8 - 3.4, 0.3, 1.5); (1.0, 2.2, 0.2 - 0.4, 1.5) 
and (1.0, 2.2, 0.3,1.2-2.3). 

Given the simplifications made above resulting in the reduction of the number of parameters from a minimum 
of 10 to only 4, the quantitative estimates obtained above for the radius, correlation length, life-time and chaoticity, 
although reasonable from the physics point of view, must not be overemphasized, the more so that the data on 
which these estimates are based, have large errors. 

^®For a more detailed discussion of this issue cf. also the comments preceding the reprinted paper by Neumeister 
et al. in ]^ ]. 
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were assumed to be real and no simultaneous fit of all correlation functions has been performed. 


Similar caveats apply to |11| where an analysis of the same NA22 data as those investigated 
in the present paper was presented and where only “marginal” agreement with the simplified 
quantum optical model of |12, was found. The fact that in the present work we could fit the 
same data supports the space-time approach and the necessity of a simultaneous fit. 

Besides the fact mentioned already that correlators in momentum-space are associated with 
a four dimensional “radius” which has no clear physical interpretation, the space-time approach 
and the momentum-space approach differ also from a purely mathematical point of view. In 
©©HI the Q dependence is postulated from the beginning while in the present paper, as in 
[p^ , it results from a complex process of integration. 

In [|^] the same momentum-space model ([^, [I^, |I^) as that used in |n| was applied and 
compared with UAl-collaboration data. This time a new technique for estimating the correlation 
data was used but only second and third order correlations were considered. At first the best fit 
parameters set for the second order correlation data was established and then used to predict the 
third order correlation function, which was found to disagree with the measured one. Given the 
insensitivity of of the fit parameters found by us on the C 2 function, this result is not surprising 
and the fitting procedures used in (as well as in in have to be qualified. Furthermore 
the general reservations expressed above about the momentum-space approach apply here again. 

One may put the question why it is necessary to invoke higher order correlations to constrain 
the source parameters whereas the presence of a gaussian density matrix implies that the first 
and second moment of the current distribution, I{k) and D{ki,k 2 ), determine all higher order 
correlations |^. 

In the applications of the space-time approach one has to take into account that due 
^"^There is of course the fact that the phases 4> of eq.(^^ enter in different combinations in the various correlation 
functions. This is true both in the momentum-space approach [l3| [l^ , as well as in the space-time approach. 
In the application to data and 0 of the quantum optical approach this complication was ’’circumvented” by 
assuming that the helds are real. 
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to the limited statistics, most experimental BEC data including those analyzed in the present 
paper give the correlations in one single variable, Q. On the other hand, the correlation function 
C 2 , for instance, depends on 6 variables: the three momenta of the pions of the pair, ki and 
k2 0. Thus a large amount of information is lost due to the integration that projects out the 
Q dependence (cf. (5) and (6)). This is a fortiori true for higher order correlations. Therefore 
phenomenologically higher order correlation data can play an important role in constraining the 
source parameters. 

On the other hand the fact that we were able to account for the data even within this 
simplified approach, proves that there are still many degrees-of-freedom not used at the present 
level of theory/data comparison so that the challenge of disproving the gaussian density matrix 
will remain a hard task for a long time ahead. 


We are indebted to W. Kittel for an instructive correspondence and to B. Schlei, U. Ornik 
and O.A.M. Helene for important discussions. N.A. was partially supported by CNPq (Brazil) 
and by FAPESP (98/01446-0). 


the case of symmetries of the source geometry, the number of independent variables is slightly reduced: 
e.g., for the boost-invariant azimuthally symmetric expanding source discussed here, there are four independent 
variables. 
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Figure Captions: 

Figure 1: The two-, three- and four-particle correlation functions in calculated using 
the space-time model defined in the text (continous lines), for parameters N2 (see Table 1) 
compared with experimental data from [pT| . 

Figure 2: The same calculation as in Figure 1 using the other set of parameters as defined 
in the text. The hgure shows that parameters which produce acceptable hts for C' 2 (Q^) do not 
necessarily give good overall hts. 


Table Captions: 

Table 1: (last three columns) for each correlation function for the parameter sets 

N1,N2 and N3. 
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Table 1 


sets 

Po 

Rchifm] 

Lr,[fm/c] 

To[fm/c] 

C2{Q^) 



N1 

1.0 

1.7 

0.2 

1.1 

1.35 

1.23 

1.49 

N2 

1.0 

2.2 

0.3 

1.5 

1.28 

1.11 

1.43 

N3 

1.0 

2.5 

0.3 

1.3 

1.31 

1.12 

1.41 
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